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February, 1994 



Progress to Date 


Substantial progress has been made toward the goals of this research effort in the past six 

months. The tasks that have been accomplished to date are: 

1. A simplified rotor model with a flexible shaft and backup bearings has been developed. The 
model is based upon the work of Ishii and Kirk. Parameter studies of the behavior of this 
model are currently being conducted. 

2. A simple rotor model which includes a flexible disk and bearings with clearance has been 
developed and the dynamics of the model investigated. The study consists of simulation work 
coupled with experimental verification. The work is documented in the attached paper. 

3. A rotor model based upon the T-501 engine has been developed which includes backup 
bearing effects. The dynamics of this model are currently being studied with the objective 
of verifying the conclusions obtained from the simpler models. 

4. Parallel simulation runs are being conducted using an ANSYS based finite element model of 
the T-501 developed by Dr. Charles Lawrence. 
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DISK/SHAFT VIBRATION INDUCED BY BEARING CLEARANCE EFFECTS: 

ANALYSIS AND EXPERIMENT 


George T. Flowers 
Fangsheng Wu 

Department of Mechanical Engineering 
Auburn University 
Auburn, AL 


Abstract 

This study presents an investigation of the 
dynamics of a rotor system with bearing clearance. 
Of particular interest is the influence of such effects 
on coupled shaft/disk vibration. Experimental 
results for a rotor system with a flexible disk 
are presented and compared to predictions from a 
simulation model. Some insights and conclusions 
are obtained with regard to the conditions under 
which such vibration may be significant. 

Nomenclature 

e — imbalance eccentricity 
f n i , f n 2 = nonlinear bearing forces 
ri — gyroscopic mass influence ratio 
i 2 = disk mass influence ratio 
xi, X 2 = shaft degrees of freedom 
x 3 , X 4 = disk degrees of freedom 
A = bearing clearance 
12 = rotor speed 

u>i = natural frequency of rotor support in Xi 
direction 

(jj 2 = natural frequency of rotor support in X 2 
direction 

u > 3 = natural frequency of rotor disk 
r = f 2 t 

£1 = damping ratio of rotor support in Xi 
direction 

= damping ratio of rotor support in X 2 
direction 

£3 = damping ratio of rotor disk 



Background and Motivation 

Rotor systems typically consist of a shaft with 
one or more bladed disks attached. Disk and blade 
vibration are issues of important concern for stress 


analysts. Complex finite element models are devel- 
oped to assess the dynamic stresses in such compo- 
nents and insure the integrity of the design. How- 
ever, models for rotordynamical analyses are typ- 
ically developed assuming that disk flexibility ef- 
fects are negligible. It is generally presumed that 
the disk is sufficiently rigid so as to not signifi- 
cantly impact rotor vibration over the operating 
region. 

There is a fairly large body of work documented 
in the literature concerning studies of disk flexi- 
bility on rotors and turbomachinery. A discussion 
of this previous work with regard to linear rotor 
systems is presented in Flowers and Ryan (1991). 
An excellent source for a comprehensive review of 
work in this area is presented by Davis (1989). 

Previous investigators into the area of coupled 
rotor/disk vibration have noted that disk flexi- 
bility has little effect on critical speeds but that 
it may significantly influence higher natural fre- 
quencies of the rotor system [Chivens and Nel- 
son (1975)]. One can draw the conclusion from 
these results that synchronous vibration due to 
imbalance will be little affected by disk flexibil- 
ity. However, there are sources of higher fre- 
quency excitation that could serve to excite these 
natural frequencies. Perhaps the most obvious 
are multi-synchronous effects corresponding to the 
blade pass frequency (from fluid forces impinging 
on the rotor blades). Another potential source is 
from nonlinearities that may be present in the sys- 
tem. For example, nonlinear bearing forces due to 
clearance effects may result in supersynchronous 
rotor vibration. There are quite a number of stud- 
ies in the literature concerning the effects of bear- 
ing clearance (and the related phenomena of rub- 
bing) on rotordynamical behavior. Some of the 
works that have most influenced the current study 
were conducted by Johnson (1962), Black (1968), 
Ehrich (1966), and Childs (1978), and Muszynska 
(1984). 



Advanced designs for many types of high speed 
rotating machinery that use magnetic bearings for 
support have been proposed and are in the de- 
velopment and construction stages. There are a 
number of such systems already in commercial 
use. Rotor systems supported by magnetic bear- 
ings must have backup bearings to provide support 
under overload conditions or if the magnetic bear- 
ing fails. Backup bearings are characterized by a 
clearance between the rotor and the bearing such 
that contact does not occur under normal operat- 
ing conditions. As a result, issues related to the 
effects of bearing clearance on rotordynamical be- 
havior are of current concern. The objective of 
the present work is to develop an understanding 
of possible coupling between the dynamics of disk 
and shaft that may be induced by bearing clear- 
ance effects and to provide guidance to designers 
concerned with such systems. 

Experimental Model 

In order to investigate whether bearing clear- 
ance can lead to significant coupling between rotor 
and disk vibration, experimental tests were per- 
formed with a rotor test rig. A drawing of the test 
rig is shown in Figure 1 . 



Horizontal Sensors 


Figure 1 Experimental Test Apparatus 

The rotor used in this study has two basic 
components: a flexible disk and a shaft. The steel 
shaft is 0.375 in. diameter and 18 inches long. The 
flexible disk is a circular aluminum plate 0.0125 
inches thick and 14.0 inches in diameter. The 
natural frequency of its lowest one nodal diameter 
bending mode is about 40 Hz. It is attached with 
bolts to a 2.0 inch diameter hub that is fixed on at 


I 


i 



Figure 2 Shaft Right End Support Device 


the center of a 0.375 inch diameter steel shaft 16.0 
inches in length. 

The rotor shaft is supported by bushings at two 
ends. The left end of the shaft is placed directly in 
the bushing base. This support provides a force 
but only a minimal couple so that the shaft is 
effectively free to pivot about this point. The right 
bearing is supported by a special device, which 
is designed to provide stiffness and to simulate 
bearing clearance, A diagram for this device is 
shown in figure 2. It has two sets of springs. 
The softer springs are to support the weight of 
the rotor, and the harder springs act as the main 
stiffness for the system. The stiffness in the x and 
y directions of each spring set may be different. 
The clearance is adjustable by turning the inner 
tapered ring in or out. In this experiment, the 
spring constant of the soft spring is 1.9 lb /in. The 
spring constant of the hard spring is 8.1 lb/in. in 
the Xi (horizontal) direction and 19.7 lb/in. in the 
X 2 (vertical) direction. 

The rotor is driven by a speed adjustable motor. 
The speed of the rotor can be controlled by turning 
the speed control knob on the front panel of 
the control box. The measurement and analysis 
system comprises proximity displacement sensors, 
sensor conditioners and a signal analyzer. A 
vertical and a horizontal displacement sensor were 
mounted to pick up the displacement signals to 
form orbit plots. A third transducer was used 






to sense the disk vibration. During the test, 
displacement signals of the shaft were sent to the 
signal analyzer, where the orbit trajectories were 
recorded and the frequency components of the 
vibration were analyzed. 

Simulation Model 

The simulation model is similar to that devel- 
oped in an earlier study [Flowers and Wu (1993)]. 
The primary differences are that rotational stiffen- 
ing of the disk have been taken into account and 
that damping is included in the bearing support 
forces. A schematic diagram of the model is shown 
in figure 3. It consists of a rigid shaft, a rigid hub, 
a flexible disk and a support with a symmetric 
clearance. The equation development is based on 
the following considerations: 

(1) The disk is assumed to flex only in the lateral 
direction. 

(2) Only rotational vibration is considered. This is 
because only the rotational motion is coupled 
with the lateral disk vibration. One nodal 
diameter disk vibration is assumed. 

(3) The rotor speed is constant. 


Flexible Disk Clearance Mechanism 



Figure 3 Simulation Model 

After some mathematical manipulations, the 
dimensionless equations of motion for the whole 
system can be obtained as 

£i + r x x 2 + r 2 x 3 + 2x 4 + f n i = cos r (1) 

x 2 - r x ±i + r 2 x A - 2 x 3 + f n2 = sinr (2) 

£3 + 2X4 + 2£3U3X3 + (u>3 — 1 )#3 

= —(£1 + 2x2) ( 3 ) 


£4 - 2x 3 + 2^3X4 + (u>l - l)x 4 

= -(£2-2x4) (4) 


fni = —<t> 


( 1 - 


fii2 = 


( 1 - 



— Ff=f )(^l X l) + 2£ 1 uqx 1 
yX] + X 2 

7 / j )(^2 x 2) + 26^2*2 

V x l + x 2 


<t>= 1 if \Jx\ +xl > S 

0 otherwise 


In the above equations, x^, x 2 , X3 and x 4 are 
scaled by e, the imbalance eccentricity. Time is 
scaled by the rotor speed, D. f nX and f n2 are the 
nonlinear moments due to contacting between the 
rotor and the bearings. 

The effective stiffness of a spinning disk can be 
strongly influenced by the spin speed. Based upon 
earlier work [Wu and Flowers (1992)], the lowest 
natural frequency of such a disk can be written as 


u 


2 

3 


19 + v ^30 
16 + a 2 


where a; 30 is the natural frequency of the nonspin- 
ning disk. This expression is used in equations 3 
and 4 to account for spin stiffening effects. 


Comparison of Experiment and Simulation 

A number of tests were conducted using the ex- 
perimental setup described in the previous section. 
In addition, studies for differing parametric con- 
figurations have been conducted using the simu- 
lation model described above. The basic results 
are very similar to those obtained by Flowers and 
Wu (1993). The results presented here are typical. 
The following discussion is directed at comparing 
the predictions of the simulation model with exper- 
imentally observed responses, with the objective of 
obtaining insight into the behavior of flexible disk 
rotor systems with bearing clearance effects. 

First, the linear characteristics of the system 
were examined. Figure 4 shows experimentally 
obtained data of the natural frequencies for the 
coupled disk/shaft vibration mode as a function 
of rotor speed. Using the data from such tests 
and additional measurements and calculations, the 
linear stiffness, damping, and mass characteristics 
of the experimental rotor system were identified. 
The numerical values are shown in Table 1. 



Parameter 

Value 

D 

0.5 

T2 

0.1 

ii 

0.06 

& 

0.06 

h 

0.01 

Vi 

0.55 


0.85 


1.8 

V 

0.3 

6 

1.2 


Table 1 Simulation Model Parameters 


Figure 4 also shows the natural frequencies 
as a function of rotor speed obtained from the 
simulation model. These results agree relatively 
well with the experimentally obtained data and 
serve to validate the structural parameters selected 
for the simulation model, ft, 3 ft, and 5ft lines are 
also shown on the figure. The intersections of these 
lines with the natural frequency curve provides 
significant insight into the nonlinear behavior of 
the flexible disk rotor system, as is discussed 
below. 



Figure 4 Rotor System Natural Frequencies 


Next, the nonlinear behavior of the system was 
investigated. A clearance value was selected and 
the rotor speed was slowly increased over a range 
of rotor speeds, with the resulting response ampli- 
tudes and frequencies measured and recorded. 

It is important to note that many vibration fre- 
quencies are possible. During the course of the 
study, lft - 6 ft frequency components were ob- 
served for certain parametric configurations. In 
addition, frequency components that are fractions 
of the rotor speed were observed for certain speed 
ranges and initial conditions. However, for a ro- 
tor with a symmetric clearance and minimal fric- 
tional effects, the odd integer multiples of the ro- 
tor speed appear to generally be the primary su- 
persynchronous frequency components. The cur- 
rent discussion will concentrate on those frequen- 
cies. However, the basic conclusions should be di- 
rectly applicable to other supersynchronous vibra- 
tion frequencies as well. 

Figure 5 shows the experimentally measured 
ampfitudes of the ft, 3ft, and 5ft components of the 
shaft response for various rotor speeds. Simulation 
results are also shown on figure 5. These results 
were obtained using the harmonic balance method 
and were verified at selected points through di- 
rect numerical integration of the equations of mo- 
tion. There is relatively good agreement between 
the simulation results and the corresponding ex- 
perimental data. Certainly the basic trends match 
quite well. 

An important observation that can be made 
from examination of these figures is that the peaks 
of the respective components correspond to the 
the intersections of the ft, 3ft, and 5ft lines with 
the natural frequency curves. It appears that the 
nonlinear effects are serving to excite the coupled 
modes of the rotor/disk system. This is true for 
both the forward and backward whirl modes of the 
rotor system. 

The peaks occurring at a rotor speed of 0.51 cor- 
respond to the intersections of the 3ft line with the 
natural frequency curve. They relate to a back- 
ward whirl mode. The peaks occurring at rotor 
speeds of 0.33 and 0.485 correspond to the intersec- 
tions of the 5ft line with the the natural frequency 
curve. Those occurring at 0.33 relate to a back- 
ward whirl mode and those occurring at 0.485 re- 
late to a forward whirl mode. Note that the ampli- 
tudes of the supersynchronous components associ- 
ated with the forward whirl mode of the disk/shaft 
vibration are much lower than those associated 
with the backward whirl mode. This result was 
observed for all the parametric configurations ex- 
amined in the course of this study. 

The ft, 3ft, and 5ft components contribute 
(more or less significantly) to the overall response 
at the corresponding intersection point. However, 
it is important to note that the ft component is 




increased rather dramatically at the intersection 
points as a result of the presence of disk flexibility. 
Apparently, the nonlinear coupling between the 
various frequency components serves to produce 
this effect. 


Simulation Model 



Figure 5. a SI Component of Shaft Vibration 
Response 


Conclusions 

A study of disk/shaft vibration induced by bear- 
ing clearance effects has been presented. Both ex- 
perimental tests and simulation results have been 
presented. The responses predicted by the sim- 
ulation model and those observed experimentally 
agree quite well. The behavior of both the ex- 
perimental test rig and the simulation model was 
quite sensitive to changes in parametric configu- 
ration. As a result it is quite difficult to predict 
exactly how a certain system is going to behave 
in the presence of a bearing clearance effect, i.e., 
whether or not coupled disk/shaft vibration will 
be significantly excited. However, a few general 
conclusions and guidelines can be drawn from this 
study. 

1. An understanding of the mechanism for cou- 
pled disk/shaft vibration induced by clearance 
effects has been obtained. The nonlinear ef- 
fects have served to produce superharmonics 
that excite coupled disk shaft modes. Con- 
versely, the additional degree of freedom pro- 


Simulation Model 



Figure 5,b 3 SI Component of Shaft Vibration 
Response 


Simulation Model 



Figure 5«c 5SI Component of Shaft Vibration 
Response 


vided by the disk flexibility has served to ex- 
aggerate the frequency components of the re- 
sponse, resulting in higher amplitudes at the ro- 
tor speeds corresponding to the respective nat- 
ural frequency /multi-synchronous line intersec- 
tion points. 

2. Bearing clearance effects primarily serve to ex- 
cite backward whirl modes for coupled disk/shaft 
vibration. Forward whirl modes may be ex- 
cited but the associated amplitudes appear to 
be much lower. 

3. The bearing support stiffnesses used in this 
study were not symmetric. These nonsymme- 
tries, together with the clearance, have served 
to produce the supersynchronous vibration that 
excites the coupled disk/shaft vibration modes. 

4. It appears to be relatively difficult to excite cou- 
pled disk/shaft vibration with bearing clearance 
effects. The occurrence of such behavior was 
very sensitive to rotor speed, as judged by the 
relative difficulty encountered during the exper- 
imental work in obtaining the peak amplitude 
responses. This conclusion is verified by the sim- 
ulation work that indicates that the behavior 
occurs only over a very limited ranges of rotor 
speed. The supersynchronous frequency must 
almost exactly coincide with a natural frequency 
in order to excite the behavior. 

5. From a design perspective, the Campbell dia- 
gram is a useful tool to predict when such be- 
havior may occur. The intersections of the cor- 
responding supersynchronous line with the cou- 
pled disk/shaft vibration mode frequency curve 
will indicate at what rotor speeds peak responses 
for that frequency component are likely to occur. 
Whether or not such behavior will actually oc- 
cur and the relative significance of such effects 
depends very strongly on the imbalance, clear- 
ance, damping, and stiffness values. 
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